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ABSTRACT 

Using prior information of population parameters based upon two auxiliary variables, a generalized family of 
efficient estimators under two-phase simple random sampling design has been proposed for estimating the population 
median of variable under study. The expressions of gain in efficiencies of proposed family of estimators with respect to 
mean squared error and bias over the existing ones have been obtained, which shows that better estimators can be 
obtained from the proposed family under the given practical situations. Numerical and graphical illustrations of the 
results obtained have also been made by taking some empirical populations considered in the literature. 
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1. INTRODUCTION 


It is well known that median is regarded as the best measure of location of population for the case of 
skewed population. In survey sampling, some authors such as Gross(1980), Kuk and Mak (1989), Singh et al 
(2001,2006), Jhajj et al (2013) etc. have defined the estimators of population median based upon variable under 
study only and using the known information of auxiliary variables highly correlated with variable under study in 
form of their parameters. In the absence of known information on parameters of auxiliary variables, two-phase 
sampling design have been widely used for developing efficient estimators of population parameters of variable 
under study practically and theoretically. 


M 7 

Assuming that information on population median z of auxiliary variable ^ (say) is known, but on 

population median ^ x of auxiliary variable A (say) is not known in a trivariate population ^ ^ which 

is estimated on the basis of first phase sample under two phase sampling design. The variable under study E (say) 
may be directly highly correlated with all auxiliary variables or may not be highly correlated with some auxiliary 
variables but correlated through other auxiliary variables. 


Considering a finite population U= (1, 2,... ^ ) in which ' and 'be the values of study 

variable E and auxiliary variables T > Z on the * unit of the population respectively. Let ^ r , ^ x and z 
be the population medians of variables E and A , Z respectively. Under two- phase sampling design, first phase 
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sample of size n is selected and variables are observed on sample units. Then second phase simple random sample of size 
m is drawn from the first phase sample and observations on the variables are noted on selected units. Let M Y , M x anc j 
M z denote the sample medians of respective variables Y , X and Y j n the second phase sample and corresponding 

M', M' ■ M' 

sample medians 7 , A and z in first phase sample. 

Suppose the values of variable Y on the sample of m units arranged in ascending order are ^ (1) ’ ( 2 )’"' , -' V ( m ) 
such that — Y(t+i) f or some integral value of t . Let P ^ m be the proportion of Y values in the sample that 

are less than or equal to the value of population median 7 (unknown). If 2 is an estimator of 1 , the sample median 
M y can be expressed as quantiles Grip) with P = 0.5. 


Defining 


p 'y = p (m x ,m y ) = 4 K (^M > with P u (x,y) = P{X < M x n 7 < M r ) 
Py- = P(m y ,m z ) = 4 { p u M - 1} > w . th P u (y, Z ) = P{Y<M Y nZ<M z ) 


P* = P, 


MJ = 4 { p « U z ) - ') ... z;, M = P( X < M, nZ < m 2 ) 


Here Pxy , Pyz and Pxz denote the correlation coefficient between estimators ^ x ^ ^ 7 ; ^ 7 ^ ^ z and 

M Y &M 7 . . 

x z respectively. 

M 

In (2006), Singh et al defined a ratio type estimator of population median 5 under two phase simple random 
sampling design using the known value of ^ z in the trivariate population ^ ^ as 


M, = M v 


M'Tf MyV( My 


M x m; 1 M z 


Where 


a t {i = 1,2,3) 


are constants. 


Up to first order of approximation, they minimized mean squared error of estimator 


mse(m s ) 


M xf{ M x) ( P^Py, ~P 
M Y f{M Y ) { p\ z -\ 
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M z f{M z )p x 


M r f{M r ) 


P xzP yz P xy 

\ Pxz 1 J 


M z f(M z )( P B P V -P^ 


M Y f(M Y ) 


V Pi ~ 1 


7 


M c 


and obtained the Bias and MSE of optimum estimator of s , given as under: 

1 


mse(m s ) 

\ / IT 


Bias(M s ) = 


7 i 

i ^ 


( 1 

1 ^ 

2 

f 1 

n 

r -»2 

— 

— 

— 

— 

Py^ 

— 

— 

R~ 

V m 

N) 


\n 

N) 

v m 

n) 

y.xz 


(1.2) 


M r 


8{M r f(M r )} 2 (l-p 2 J 2 

K 


1 1 w 2 

' Pxy ~ Pyz Pxz) ~ 2 Pxy ( P xy ~PyzPxz ){ 1 ~ Pxz) 


/ \ / 2 \| 
M x f(M x y Pxy PyzPia P 

~ y J j ( P yi — Pxy Pxz )” ~^^Pxz ( P xy — P yzP . xz ) ( Pyz ~ P xyPxz} 


1 1 
m 


(1.3) 


R 1 = 

y.xz 


Where 


Plv + Pyz ~ ^P xyPyzPx 

!- pi 


Using the knowledge on range of variable Z in addition to its population median ^ z , Gupta et al (2008) 

M 

defined the estimator of population median Y under the same sampling design considered by Singh et al (2006) as 


M p = My 


[ftp 

r i 

[m z +rP 

Yi 

[m z + rA 



V M Z + R Z y 


v Mz 


Where 


Yt (i = 1A3) 


are constants. 


(1.4) 


M, 


Up to the first order of approximation, they minimized MSE of p for 


Y i = - 





1 

Mz f(M z )p„ 

Pxz Pyz Pxy 

M x f(M x ) 

( P xzP yz Pxy 

/ 2 

f M? ] 

Myf{My) 

\ Pxz 1 y 

M Y f{M Y ) 

^ Pxz 1 y 


y M Z + R Z y 
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n=- 


M z f{M z ) 


f PxzPxy~Pyz^ 


M, 


V J 


M Y f{M Y ) 


\ y.xz 


pi - 1 


M 

And obtained the bias and mse of optimum estimator of p , given as under : 

M v 


Bias(M p ) = 


8 {M Y f(M Y )} 2 (1 - pl z ) 2 

X ( m ~ n ) { ( Pxy ~ PyzPxz ^ “ 2pxy ( Pxy “ PyzPxz X 1 ~ ) 

M y f (M y) / W 2 \ 1 

M x f(M x )^ Pxy P >-- P ^ Px *)\ 

+ ( m ~ I ( Pyz ~ PxyPxz X + 2Pxz ( P ' y ~ PyzPxz ^ Pyz ~ PxyPsz ) 

{pyz~ PxyPxz)^ pj)} 


Myf(M y ) 


m 7 


V M z + Rz 


M z f{M z )\ 

+ ( “ "^){ P * yP., P::P.S ~PP,P,: {p V ~ PyJ>^~ Pi) 

p.{p, /\/\)c pi)} 


Myf{M r ) 


M z f (M z ) 


' M z A 
yM z +R Z J 


(m p ) - 1 

Yi-i) 


12 f 1 0 

p.._ — 

\r* 

v 4{f(M Y )} 

Y m N ) 

U N ) 

' yz 

1 ' 1^777 77 j 



mse(m s ) 


(1.5) 


(1.6) 


Recently Jhajj et al (2014) defined an efficient family of estimators of median using the same information used by 
Singh et al (2006) under same sampling design as 


M m =M y 


f M- A 




f M. V2 


\ M U 


'3 \ 1VI X 1V1 X 

M y +M y 


(1.7) 


Here ^ >1 , ^ 2 and are constants. 


M . 

Up to the first order of approximation, they minimized MSE of YH for 




Myf(M v ) 


PxzPyz P xy 


V Fxz 


pi - 1 


M z f{M z ) 


M Y f{M Y ) 


PxyPx: Pyz 


V P’xz 


pi - 1 
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2 M x f(M x ) 


M Y f{M Y ) 


PvzPxz - Pxy 

Pi~ 1 


V 


y 


rM v 


And obtained the bias and MSE of optimum estimator of YH , given as under: 
M v 


Bias(M iH ) = 


8 {M r f (M t )} 2 (1 — pl z ) 2 

: im' “ n ) { ~ P * P = ^ ” 2px * “ P ’- P - ) (* " p i ) 

+ l ( P * ~ P " P = ) 2 + 2p - ( P * “ p p ' H/’,' - PxyPxz ) 

- 2 p y: [p r _ - p xfPxz )(l ~pl)+ ^j^ {Pr- ~ P*yP*-- ) (' - P~ ) } 

+ ( ~ n - { P= K ~ PyJ>Jf - 2P^, {p* - /’ /’■)(' Pi) 


( M \ ~ 1 

Y i n 


fi l > 

Pi~ 

f 1 

n 

K. 

1 YH min ~ 4{/(M y )} 2 

y m N ) 


,n N j 

K m 

nj 


( 1 . 8 ) 


(1.9) 

mse(m s ) 

2. PROPOSED ESTIMATOR AND ITS RESULTS 

Suppose that variable X is highly correlated with variables Y arK j Z jjj the trivariate distribution of variables 
Y , X anc j Z > but variables Y an( j Z are correlated to each other only because of their correlation with X when 

population median ^ z of variable Z ; s 0 nly known, we propose a new generalized family of estimators of population 

M 

median ' under two phase sampling design described in Section 1 as 


^ me ~ 




f 2rf\M x -M' x )\ 


m 7 

72 

m 7 

m y +g(m'-m y 

exp 






L V / J 


{ (m x +m' x ) J 


l M 'z\ 


V M 'z\ 


( 2 . 1 ) 


Where are constants with P > ® . 
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M 

Assuming that mean squared error of the estimator me exists and finite. So, for finding its mean squared error 

and bias up to first order of approximation, we expand j n terms of s s . up to their second degree by defining £ s 

as 


My . My , My - 

s n = — - - 1 s, = — - - 1 e 7 = — — - 1 

My My My 


M 'x 1 

£, = — — - 1 

M v 


M z M' 

£, = — - - 1 £,= — - - 1 

4 M z M z 


provided ^ , such that ^( £ i) ^ f or * = 0, 1, 2, 3, 4, 5 


M 

Up to first order of approximation, bias and mean squared error of iHe obtained are : 


Bias(M m0 )=— 


flj} 

ytl TV 




Mm z) ) 




\m n 

( 

+M 


)\ 2 v " v "" 2 

(M x f{M x )Y {M z j\M z )) 


-A 


m 


(Myf{M T )) +ThP » (Myf(My)) 


+ mA 


(M x f{M x ))~ 


(M z f{M z )) 


(2.2) 


MSE 


(^me\ 


Ml 


1 1 


--- (M Y f(My)] 


(\ n 


~ 2r hpy 


(M z f{M z )Y 


+ 


\ n TV j 

{--- 

v m n 


{Myf(My)) 

+rj{ (M x f(M x )f +v 2 (M z f(M z )y +2^ (1 -Q)p 


ril(M z f(M z )) 

-^\{(e 2 -7B)(M Y f(M Y )) 2 

{M x f(M x )j' 


(M r f{M r )) 




(M r f(My)) 


(M z f(M z )) 


a MSE[My m ). 

ror any nxed value ot , x 7 i 


^i = (l-6>) 


PyzPxz P xy 


1 -P 


is minimized for 

M x f(M x ) 


1 


M r f(M r ) 


(2.3) 
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( 'N 

r xyrxz r yz 

mj(m 2 ) 

V * ~ Pxz y 

M v f(M v ) 


, M 2 f(M 2 ) 
V M r f(M r ) 


And its minimum value is given by 


MSE(M me )-- 


1 


“ ¥ H K -1 III 1 -")’ Ol.-f - 2 ") 


m IV 1 1 77 N 


1 1 


4 {f{M T )) 

And the bias of its optimum estimator obtained is 

1 


1 1 


Bias(M me ) — 


[M 

1 

1 

\m n ^ 

V L 


+ 2 P„(p,,-Az -Pxy)( 1 -P^)+ 2 P^(P^ -P,,)(!-P=) 

+ 2 P„ (P„P C -Pvz)(p„Pv_- -Pr,)+(P.,.P= -Pv--)* 

-(p P -P )( lV) ^/K) ( ) (lV ) ^/W 


+ \ X -hr, 


M x f{M x ) 
M Y f(M r ) 


n Nj^{M z f(M z ) 


~P, 


(2.4) 


(2-5) 


3. COMPARISON 


Since the mean squared errors of optimum estimators defined by Singh et al (2006), Gupta et al (2008) and Jhajj 
et al (2014) is same up to first order of approximation, so we compare the proposed family of estimators with Singh et al 
(2006) estimators. 

Using the expressions (1.2) and (2.4), we have 


MSE(M s ) mn -MSE(M me ) 


1 

4 (f(M Y )f 


m 





> 0 for 0 < 6 < 2 ^ ^ 

From (3.1), we note that the proposed family of estimators is always better than Singh et al (2006) estimator, 
Gupta et al (2008) estimator and Jhajj et al estimator (2014) for ^ < 9 < 2 

From the expressions of the biases obtained, we noted that no concrete interpretation can be given theoretically 
about the superiority of one estimator over the other on the basis of their biases because of the complicated expressions 
obtained in comparisons. To have a rough idea, we compare the biases numerically and graphically by taking data of 
populations considered in literature in Section 4. 
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4. NUMERICAL ILLUSTRATION 

For illustration of results numerically and graphically, we take the following data 

Data 1 (Source: Singh, 2003). Y: the number of fish caught by marine recreational fishermen in 1995; X: 
The number of fish caught by marine recreational fishermen in 1994; Z: The number of fish caught by marine recreational 
fishermen in 1993. 


N = 69, 

^=0.1505, 

My =2068, 

} ') = 0.00014 

n = 24, 

Pyz =0.3166, 

M *=2011, 

/ ( M ^) = 0.00014 

m =17, 

Pxz = 0.1431, 

Mz =2307, 

/( M z) = 0.00013 


Data 2 (Source: Aezel and Sounderpandian, 2004). Y: The U.S. exports to Singapore in billions of Singapore 
dollars; X: The money supply figures in billions of Singapore dollars; Z: The local supply in U.S. dollors. 


so 

II 

£ 

Pxy = 0.6624, 

=4.8, 0.0763 

n = 23, 

Pyz = 0.8624, 

Mx =7.0, f ^ Mx ) = 0.0526 

m =15, 

Pxz = 0.7592, 

^ z =151, f( Mz ) = 0.0024 


Data 3 (Source: MFA, 2004). Y: District-wise tomato production (tones) in 2003; X: District-wise tomato 
production (tones) in 2002; Z: District-wise tomato production (tones) in 2001. 


N = 97, 

Pxy = 0.2096, 

My =1242, 

f y) = 0 00021 

n = 46, 

Pyz = 0.1233, 

Mx =1233, 

f{ M x) = 0000 22 

m =33, 

Pxz = 0.1496, 

Mz =mi , 

f { M z) = 00002 3 

For Data 1 





Table 4.1: Bias and Relative Efficiency of Estimators 




Bias 

Relative Efficiency 

M s 

M P 

M m 

M YH e 

M s 

M me 

-0.13 

58.62 

20.22 

32.57 

16.28029 

100 

90.34 

-0.08 

58.62 

20.22 

32.57 

16.64022 

100 

93.99 

-0.03 

58.62 

20.22 

32.57 

16.97063 

100 

97.72 

0 

58.62 

20.22 

32.57 

17.15472 

100 

100 

0.3 

58.62 

20.22 

32.57 

18.41137 

100 

124.36 

0.6 

58.62 

20.22 

32.57 

18.60579 

100 

147.62 

0.9 

58.62 

20.22 

32.57 

17.73798 

100 

161.34 

1.2 

58.62 

20.22 

32.57 

15.80795 

100 

158.39 
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Table 4.1: Contd., 

1.5 

58.62 

20.22 

32.57 

12.8157 

100 

140.45 

1.8 

58.62 

20.22 

32.57 

8.76122 

100 

116.04 

2.1 

58.62 

20.22 

32.57 

3.644517 

100 

92.54 


For Data 2 




0 


Figure 4.2 


Table 4.2: Bias and Relative Efficiency of Estimators 


e 

Bias 

Relative Efficiency 

M s 

M p 

M m 

M Y hb 

M s 

M Y H0 

-1.7 

0.37 

0.28 

0.026 

0.30093 

100 

26.19112 

-1.2 

0.37 

0.28 

0.026 

0.15731 

100 

36.759 

-0.7 

0.37 

0.28 

0.026 

0.05226 

100 

54.14859 

-0.2 

0.37 

0.28 

0.026 

0.014208 

100 

83.533 

0 

0.37 

0.28 

0.026 

0.029993 

100 

100 
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Table 4.2: Contd., 

0.4 

0.37 

0.28 

0.026 

0.043043 

100 

140.2007 

0.8 

0.37 

0.28 

0.026 

0.031403 

100 

175.4711 

1.2 

0.37 

0.28 

0.026 

0.00493 

100 

175.4711 

1.6 

0.37 

0.28 

0.026 

0.06595 

100 

140.2007 

2 

0.37 

0.28 

0.026 

0.15166 

100 

100 

2.4 

0.37 

0.28 

0.026 

0.26207 

100 

69.92491 

2.8 

0.37 

0.28 

0.026 

0.39716 

100 

49.91067 

3.2 

0.37 

0.28 

0.026 

0.55695 

100 

36.759 

3.6 

0.37 

0.28 

0.026 

0.30093 

100 

27.92802 

4 

0.37 

0.28 

0.026 

0.15731 

100 

21.814 



Figure 4.3 



Figure 4.4 
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Table 4.3: Bias and Relative Efficiency of Estimators 




Figure 4.5 



0 

Figure 4.6 
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From tables (4.1), (4.2) and (4.3), we noted that proposed family of estimators have smaller mean squared error 
than Singh et al (2006) estimator, Gupta et al (2008) estimator and Jhajj et al estimator (2014) as well as have smaller bias 
than Singh et al (2006) estimator, Gupta et al (2008) estimator and Jhajj et al estimator (2014) in all the three populations 

<9e(0,2) 

considered corresponding to v ' . The graphical representation also supports the numerical results. 

CONCLUSIONS 


1/ Cl V7, Z- I . . 

It is interesting to note that for any value of v ' , the proposed sampling strategy for estimating the 
population median is efficient than the sampling strategies obtained by Singh et al (2006) estimator, Gupta et al (2008) 
estimator and Jhajj et al estimator (2014). Numerical results obtained corresponding to all populations considered in 
literature also justify that proposed sampling strategy of estimating population median having smaller bias is efficient than 
existing ones. Comparison of the estimators is shown graphically based upon the numerical results obtained also give same 
type of trend. Therefore, we strongly recommended that proposed sampling strategy should be used for getting the accurate 
value of population median. 
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